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Abstract 

Our recent method to calculate renormalized functional determinants, the partial wave cutoff 
method, is extended for the evaluation of 4-D fermion one-loop effective action with arbitrary mass 
in certain types of radially symmetric, non-Abelian, background gauge fields (including instanton- 
like and instanton-antiinstanton-like configurations). A detailed study on functional determinants 
for matrix-valued radial differential operators is presented, explicating both our analytic treatment 
on the high partial wave contribution and the application of the generalized Gel'fand-Yaglom 
formula to determine the low partial wave contribution. In general, some numerical work is needed 
for the low partial wave part. In the massless limit, however, the factorizable nature of our partial- 
wave radial differential operators can be exploited to evaluate semi-analytically even the low partial 
wave part, and we thus have the full fermion effective action calculated explicitly in a class of non- 
Abelian background gauge fields. With nonzero mass, we also perform necessary numerical analysis 
as regards the low partial wave contribution to produce numerically exact results for the massive 
effective action. Comparing these against the results of the large mass expansion, the validity 
range of the large mass expansion is addressed. Also studied is the fermion mass dependence of 
the effective instanton-antiinstanton interaction. 
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I. INTRODUCTION 



While in field theoretic studies we are often led to consider the one-loop effective action 
in some nontrivial backgrounds, it is quite difficult to have it explicitly evaluated. Also 
lacking are well-controlled approximation schemes for the quantity, which can cover broad 
types of backgrounds. This is true in four spacetime dimensions especially. Recently, for the 
case involving radially symmetric backgrounds, we (with G. Dunne) [U [2] developed a new 
partial-wave-based calculational scheme, the partial wave cutoff method, by which the exact 
computation of fully renormalized one-loop effective actions can be performed explicitly. 
This method is a unique package of analytical and numerical procedures (to treat high and 
low partial- wave contributions, respectively). So far it has been applied to the accurate 
determination of QCD single-instanton determinants for arbitrary quark mass values [3] , to 
the prefactor calculation in the false vacuum decay [3], and to the evaluation of the scalar 
one-loop effective actions (for any given mass values) in a class of Abelian or non-Abelian 
radially symmetric background gauge fields [2]. Also, very recently, the fermion one-loop 
effective action in Abelian radial background gauge fields has been studied by this method 
[5J, an important byproduct of this work being that there exist marked differences between 
the small mass mass limits of the derivative expansion for spinor and scalar theories. 

In this paper the partial wave cutoff method will be used to study 4-D fermion one- 
loop effective actions in a class of genuinely non-Abelian, radially symmetric, background 
gauge fields. This case differs from those of our earlier studies in that, as the differential 
operators pertaining to partial wave sectors are not completely separate, we here have to 
deal with an infinite number of functional determinants for matrix-valued radial differential 
operators. One might then suspect that, because of technical difficulties in renormalizing 
the infinite product of such functional determinants and also in performing the needed 
numerical calculations, our whole approach becomes impractical in this case. Despite this 
complication, it will be demonstrated here that our method can suitably be extended such 
that the exact computation of the effective action becomes possible for this well. 
There are also issues specifically involving fermion effective action (e.g., the massless limit 
behavior), and we intend to provide clarification on such aspect. 

Specifically, in a 4-D Euclidean SU(2) gauge theory, we will consider in this work the one- 
loop effective action of a Dirac field (with mass m and in the fundamental representation) 
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when the background gauge fields are given as 

A,(x)=r ] ^ a x u f(r)r a , (1.1) 

where p, u — 1, 2, 3, 4, r = ^Jx^x^ and r/W are 't Hooft symbols |6j. Writing /(r) = ^H(r), 
the radial function H(r) is then assumed to have the form 

(Case I):g I (r;a)= 1 j^fL , (H > 1) , (1-2) 

f r /n)2 l + tanh(^) 

(Case II) : # n (r; i2, /3) = 1 '/ j .„ f/3 can take either sign) (1.3) 

l + ( r /P) 2 

with constant background parameters p, a, R and 0. This genuinely non-Abelian back- 
ground field has been chosen so that one can learn something about the behavior of the cor- 
responding fermion effective action as one changes the background parameters (and fermion 
mass m). Note that, in Case I (and Case II with (3 > 0), our background gauge fields have 
the Pontryagin index equal to ±1; on the other hand, Case II with /3 < (and the ratio R/p 
significantly larger than 1) corresponds to a well-separated instanton-antiinstanton configu- 



ration with zero Pontryagin index. [The Pontryagin index of the fields (1.1) is determined 



by the two numbers (see ( 2.13[ )), H(0) and H(oo)]. In Fig. [T]we have given the plots of the 



function H(r) for some representative choices of our free parameters. Needless to say, with 
these backgrounds, the small mass limit of the fermion one-loop effective action becomes 
particularly interesting physically (because of the issue concerning fermion zero modes). 
In a well-separated instanton-antiinstanton configuration, the related effect is also believed 
to generate long-range interaction between instanton and antiinstanton [7J [8]. Our study 
should illuminate such aspect, too. 



For the choice a = +1 or —1 in Case I above, our fields (1.1) represent single instanton 
or antiinstanton solutions [9], in the regular (for a = +1) or singular (for a = —1) gauge. 
These are (anti-)self-dual backgrounds, and here we have a simple relationship [10J between 
the fermion one-loop effective action and that for a scalar field, which can be exploited for the 
effective action calculation. This was crucial in the calculation of [3]. But, for other cases we 
will consider, we need to perform an entirely separate calculation for fermions; in this direct 
fermion analysis, additional complication due to the magnetic moment coupling term arises. 
Here we have found that the recently established chiral separation of the fermion effective 
action jTT] can be utilized to the full advantage — thanks to the latter, our technique to 




(c) (d) 

FIG. 1. Plots of the radial profile function H(r). For Case I, we have drawn the profile for 
a = 1,3,5 in (a) and a = —1,-3,-5 in (b). Note that H(r) behaves like a step function if \a\ 
becomes very large. The plot in (c) — a spherical-wall-like (anti-)instanton configuration — is 
appropriate to Case II with /3 > 0, and the plot in (d) — an instanton-antiinstanton configuration 
— corresponds to case II with (3 < 0. 



evaluate functional determinants of radially separable differential operators can be extended 
to our problem involving Dirac fields. 

This paper is organized as follows. In Sec. [IT] we describe a general outline for the cal- 
culational scheme we shall use, and also collect, for later use in the paper, various useful 



formulas. This is followed in Sec. Ill by our analysis for the high partial-wave contribution 
to the fermion effective action; this part is calculated analytically using a WKB-type asymp- 
totic series, for matrix-valued radial differential operators in this paper. We then combine 



this high partial-wave contribution with the low partial-wave contribution in Sec. IV, to 



obtain the explicit results for the renormalized fermion effective action with the background 
fields specified as above. An important part in the computation of the low partial wave 
contribution is the evaluation of functional determinants with matrix- valued differential op- 
erators, and for this we use the (generalized) Gel'fand-Yaglom approach [12] . The partial 
wave contributions must be summed according to certain specific grouping procedure. In 
the limit the fermion mass approaches zero, we can give a semi-analytical treatment for 
partial-wave functional determinants and the results are used to discuss various features 
that the fermion effective action in the massless limit exhibits. (See for instance our formula 



(4.58), where the exact small-mass-limit form of the fermion one-loop effective action in our 
general Case I backgrounds is given). This is possible because, unlike the scalar field case, 
partial- wave radial differential operators in the spinor case enjoys certain factorization prop- 
erty. Also presented are the results with nonzero mass. This involves an extensive numerical 
work, but, with our method of accelerating the convergence, numerically exact effective ac- 
tion can be obtained without excessive labor. For large fermion mass, we recover from our 
expression the result of large mass expansion. From our computation of the fermion effective 
action in Case II backgrounds, certain information on the effective instanton-antiinstanton 
interaction can also be gained. Sec. M contains concluding remarks. 



Some supplementary results, related to the high partial-wave contribution of Sec. Ill 



can be found in Appendix [A] In Appendix [B] we clarify some subtle aspect arising when one 
uses the Gel'fand-Yaglom approach in computing partial-wave functional determinants in 
the massless limit. 



II. PREPARATORY SETUP FOR OUR COMPUTATION 
A. Fermion effective action in radial backgrounds 

The bare fermion effective action is 



T(A] m) ~ — lndetf— ■ D + m] + const. 

~ — lndet[(7 ■ D) 2 + m 2 } + const., 
2 
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(2.1) 



where 7 • D = 7 M -D M , Z)^ = — iA^, and {7^,7^} = — 25 Ml/ . Its Pauli-Villars regularized 
form, using the Schwinger proper-time representation [14] , is given by 

1 f°° ds / 



F(s) 



-A 2 1 



F s) 



d x Tr ( x 



-sin-D) 2 



-s -d 2 



x 



(2.2) 
(2.3) 



where 'Tr' denotes the trace over Dirac spinor and internal isospin indices. Then the renor- 
malized fermion effective action in the 'minimal' subtraction scheme can be identified with 
the expression 

I I . A 2 

J 2 



T zen (A; m) = lim { T A (A; m, 



A 2 r 

In — J d A x tr(F^F ; 



\W ) < 1 



(2.4) 



where /1 is the renormalization scale, 'tr' denotes the trace over isospin indices only, and 
F^v = i[D/j,, D u ] = d,j,A v — d u A^ — i[A^, A„]. For large mass m the DeWitt WKB expansion 
(or heat kernel expansion) can be used to generate the large-mass approximate form for the 
effective action. But, in this paper, we are more interested in the exact evaluation of the 



effective action in given backgrounds of the form (1.1) (but mass kept to arbitrary value). 



To compute the effective action with the 'radial' background (1.1 ), it is convenient to use 
the chiral representation for 7-matrices 
/ 



7m 

We then have 



V 



a„ 







(with cr M = (a, %) and 07, = (a, -i) = (crD) 



(2.5) 



(l-D) 



_n2 _ 1 (-) F 






_n2 _ l rj (+) . p 



(2.6) 



2 'llivaya,*- fiu 

(here D 2 = D^D^), and as a result the fermion effective action can be expressed by the sum 
of chirally projected ones pj], viz., 

r rcn (A; m) = rW(>l; m) + r£(A; m) (2.7) 

with 

r°° ds 



m) 



lim - 

A^oo 2 



(V™ 2 * - e ~ A2s ) F^\s) 



2 s -A 2 s 

- [ c — e 

s 

1 1 In ^ J d 4 x (tr {F^F^) T ~tr {F^F, 



3(4tt) 



F^(s) = J 



d x Tr ( x 



-D 2 — rr F 



-s(-d 2 ) 



X 



(2.8) 
(2.9) 



where *F fJiU = \e^ v xsFxs and 'Tr' denotes the trace over 2-component (i.e., in the given chiral 
sector) spinor indices and isospin indices. Note that, to renormalize T^, we need to include 
the term proportional to the Pontryagin index in addition to the term involving the classical 
Yang-Mills action. To obtain the full effective action, there is no need to compute and 
T^~2 separately, the two being related by pj] 

r<+>(4 m) - m) = In ^ J d 4 x tr(F, u *F, u ). (2.10) 

Explicit evaluation for one, the simpler from the two quantities fW and for a given back- 



ground field, thus suffices. For our background fields (1.1), i.e., for A^(x) = viiil x vf(r)r a , a 



simple calculation shows that r)ffi a o - a F \ v (but not r]^) a a a F^ u ) takes a purely radial form 



n^rr F 



fiua 
2 rf„\2 



4/(r) + rf'{r) -2r 2 f{r) 



<J a T a 



-2g F (r)a a T a , 



(2.11) 



i.e., a radially separable differential operator is obtained only for a particular chiral com- 
ponent of the fermion quadratic operator (7 • D) 2 . Hence, for our field A^(x) with the 
^(+)_ S y m bol chosen, it is the quantity F^~2 that we may try to evaluate directly by applying 



our partial wave cutoff method; for PW, on the other hand, we can use (2.10). So the full 



effective action r ren follows from the result for r{~) alone. The situation is just the opposite 
if our field A^x) happens to involve the f/~)-symbol. [This choice of depending on our 
background field form is appropriate with an arbitrary radial function f(r); for a particular 
form of f(r) which gives rise to (anti-) self-dual field strengths, the other choice would be 
more suitable (with rj^laaF^ = for the corresponding backgrounds)]. 



When the background field is given by the form (1.1), we have the classical action ex- 
pressed using the function H(r) = r 2 f(r) as 

dr 

/} I r/'.r tr/- - I ' / 

and the Pontryagin index as 

1 

16tt 2 



l - J d A x trF 2 u = 12tt 2 J™ y {r 2 H'{rf + AH (r) 2 [H(r) - l] 2 } , (2.12) 



w 



*F 



2H(rY - 3H(r) 



r=oo 
r=0 



(2.13) 



Hence, if the ^"^-symbol is chosen in our expression for A^(x), we find for the two Cases in 



(1.2) and (1.3) the Pontryagin index 

(Case I) : w 



1 

-1 



a > 1 
a < -1 



(2.14) 



(Case II) : w 




(2.15) 



Also, for the background field form (1.1 ), the differential operator we must deal with, — D 2 
l^u^aCa-^ty; can be expressed in the form 



-D 



n^cr F 



d 2 



dr 2 



1^. + ^ + 8 /( r ) f • li±) + 3r V(r) 2 + Ag F {r)S ■ f, 



where T a = |r a , S a = ^a a , and and L 2 are specified as [6] 



£,(±) = -l„ (±) x d 



r (±) r(±) 



l^abc 



(2.16) 



(2.17) 



For the evaluation of rW(v4;m), we may then resort to a kind of block diagonalization for 



the differential operator (2.16) in the form of partial waves. See Part B of this section.. 



Another useful information as regards our background field form ( 1.1 ) is that the appear 



ance of the ?p + )- or rf ^-symbol in the expression is actually tied up with the gauge choice. 
Explicitly, using the relation 



(±)^a 



;T a = iQ { ^(x)d^Q i±) (x), 
it is not difficult to show that 



n 



(±)W 



24 =F ix ■ r 



r 



(2.18) 



n 



(±)W 



I L Xy 



#(r)r a 



^(x) + in (±) (x)a M n^(x) = ijgj^i - #(r)]r B , (2.19) 



i.e., under the (singular) gauge transformation involving the SU(2) matrix f^Wrr), our 
^w-symbol form with the radial function H(r) goes over to the r/^-symbol form with 
the radial function 1 — H(r). This has the consequence that an identical fermion effective 
action T ren (A; m) would result either with the form A^ix) = rj^ a ^H{r)T a or with the form 
A^x) = i]j^l%H(r)T a where H{r) = 1 — H(r). Then notice that, for H(r) given by the 



form in (1.2), we find 



Hi(r) = 1 - 



(r/p) 



2a 



{r/p) 



-2a 



(2.20) 



l + (r/p) 2a l + (r/p)- 2a ' 
Therefore, in our Case I we may just evaluate the effective action with the radial function 
of the background field chosen as 

(r/p) 2 " 



Hi(r;a) 



(r/p) 



2a ' 



(«>i); 



(2.21) 



1 | — 



r/p 



Rip 



r/p 



(a) 



(b) 



FIG. 2. Plots of the radial profile function Hu(r) given in (2.22): (a) for (3 > and (b) for (3 < 0. 



then the gauge invariance of the effective action tells us the result appropriate to the same 
function form for H(r) but with a < —1. Also, applying the above gauge transformation to 



our Case II, we may replace our radial profile function in (j 1 . 3p by the form 

(r/p) 2 



H n (r;R,P) = l- 



1 + (r/p) 2 
r/p)- 2 



(r/pY 
1 + (r/pf 

+ H n (r;R,-(3), 



1 



tanh( r - K 



0P 



(2.22) 



1 + (r/p)- 2 

as this form should lead to the same set of effective actions. Note that for (5 < and R large, 



what we have in (2.22) is the sum of a singular-gauge instanton (antiinstanton) located near 



the origin and a spherical- wall-like antiinstanton (instanton) configuration at large radius R. 



See Fig 2 for the illustration of the radial profile function corresponding to our form (2.22) 



B. Partial- wave decomposed form 



The differential operator in (2.16) can be decomposed into an infinite number of partial- 



wave radial differential operators (with matrix coefficients). For our partial waves, let us 
consider the basis \j, j 3 , q, I, l 3 ) where various quantum numbers introduced are specified by 



(L 2 )' = 1(1 + 1), 
(Q 2 )' = q(q+l), 

(j*y=j(j+i), 



1 3 

/ = 0, -, 1, -, 
'2 2 



(with Q a = L^+T a ), q 
(with J a = Q a + S a ), j = 



l± 



(2.23) 



W = ja = ~3, -3 + 1, 



(L 



(±)y 



Z 3 = -/,-/ + !,-••,/. 



In this basis, the operator in (2.16) is not completely diagonal, but we may still write it (for 
given values of / and j) as 



D 



?7 (±) er F 



H 



d 2 3d 4/(Z + 1) A r . . 



+3r 2 f(r) 2 + 4g F (r)S -f 



= -% + V y (r) 



g(g + !)-/(/ + !) 



(2.24) 



with suitable matrix 5 • T in the space of allowed g-states (for given I, j). [In (2.24), 
•9( 2 «) = ^5" + f — ~2 represents the 4-D Laplacian c^c^ for given angular momentum]. 
As for the matrix S ■ T, we here note that, for given I, the quantum number q should be 
equal to I + - when j = Z + 1 and equal to / — | when j — I — 1, while g can take either value 
of Z ± 2 when j = I ^ 0. When j = / = 0, only the value q = | is available. Then, after 
somewhat lengthy but straightforward calculations, we obtain following representations for 
S-f: 



S-T 

S-f 
S-f 



1 

4 



, if j = I ± 1 ^and q = I ± - 



1 



-2/ -3 4 v / /(/ + l) 



if j = /^0, 



(2.25) 



if j = / = ( and g 



1 



4(2/ + 1) [ A^Ul + T) -2/ + 1 
3 

where our 2x2 matrix form for j = I ^ is written relative to the basis (|g = / + ~), |g 



/ — g)). Using the expression in (2.25), we thus find, for the 'potential' V/j(r) defined in 



(2.24) 



V lH1 (r) = 3r 2 f(r) 2 + Alf(r)+g F (r), 
Vy_i(r) = 3r 2 /(r) 2 - 4(Z + l)/(r) + ^(r), 



(2.26) 
(2.27) 



V M (r) = 3r 2 /(r) 2 + 4/(r) 



/ 



/ 







-/ - 1 



+ 



^(r) -2/ -3 4 v //(/ + l) 



2/ 



1 V 4 V /Z(/ + 1) -2/ + 1 
(Z^O) 



with # F (r) = 4/(r) + r/'(r) - 2r 2 /(r) 2 (see pll] )), and 
V ,o = 3r 2 /(r) 2 -3( ?F (r). 



(2.28) 



(2.29) 
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[The bold-faced letter for Vi,i(r) is to indicate that it is matrix- valued] . To each partial wave 
labeled by quantum numbers J = (l,j,j3,h) corresponds the (matrix) radial differential 



operator Uu specified by (2.24) and (2.26)-(2.29). 



Our formula (2.8) for the effective action T^(A;m) can then be reexpressed using the 
quantity involving the radial differential operator Tiij appropriate to partial waves. Here 
note that, as should be evident from the discussion above, fW(i;m) in the background 



(1.1 ) with the rf ^-symbol picked has the same value as T^(A; m) in the very background 



but with the 7/ + )-symbol taken. Then, thanks to (2.10) and (2.13), we are led to conclude 



that the full effective action T ien (A;m) in the background (1.1) is represented by the same 
function of f(r) irrespectively of which //-symbol is picked. Knowing this, we may well 
consider only the background field form 



(2.30) 



from now on and go on to evaluate the quantity r[ e 2(A;m), a particular chiral projection 



of T Ien (A;m), in this background. The radial function H(r) is that of the form (2.21) or of 
( 1.3[ ). Note that, because of (2.10), the full renormalized effective action can then be found 
simply by using the formula 

r ren (A; m) = 2lt n )(A; m) + ^ / d 4 x tr(F^F^). (2.31) 

Now, for the partial- wave-decomposed form of r^~2(A;m), we may express the function 



F^(s) (see Q) by the form 

j 

F&\s) = J^dr tr{Gy(r,r; a ) - Gf ee (r, r; S )} 
where we introduced the radial proper-time Green's functions 



sH 



G\ 



Gi,j(r,r';s) = 
of the operators "Hy and ois defined through [I] 

^■1,3 = ^ n i,i r3/2 



free , 



r, r ; s 



(2.32) 



(2.33) 



£ 41(1 + 1) + * , 
~d^ + * + V ^ (r) ( = 



(2.34) 
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When I is large, we can employ the WKB method to calculate f[j \s) systematically; 



this, we do in Sec. Ill (Using the result, the renormalization problem is also solved in an 
expedient way). But we need a different strategy for small I, and to determine the related 
contribution to T(. a }(A;m) it is more convenient to consider the quantity resulting after 
s-integration, 



f 



ds 



In 



det(Kzj + m z 
det(— Sns + to 



In 



det(Hi tj + m 2 ) 
det(— d 2 ^ + to 2 ) 



(2.35) 



Then, given a specific background, we may evaluate the latter quantity with the help of the 



Gel'fand-Yaglom method. We do this for our (matrix) differential operators Tiij in Sec. IV 



In performing these calculations, some care must be exercised if the effective radial potential 

V l , J (r) = m + l ) + 3 ' +V l Ar) (2.36) 



doest not have the same small-r and large-r behavior as V l 



free / 



«l±3/4. i n the latter 



case, to avoid the appearance of largely oscillating terms in the process of summing partial 
wave contributions, a certain grouping of terms might be contemplated [5]; this is relevant 
also with our background fields since, or our Case I for instance, we have /(r) ^ 



and so 



Vu(r) ~ 4 ^^ 2 )+ 3 / 4 as r -> oo. 



To facilitate our discussions, we will now write out the mathematical expressions we must 
evaluate to determine T^~2(A; to) for our background fields. Following Refs. [3], we introduce 
the (floating) parameter L as the partial wave cutoff and write 



r( - )f a- 

1 ren 



lim 

L— ¥oo 



r\-l(A;m)+r\;> L (A;m 



(-) 



(2.37) 



here, T^^A; to) represents the 'low' partial wave contribution given as 



KL 



(A;m) 




det("H ,o + rn? 
det(-<9 ( 



(1=0) 



+ m z 



+ E |(2i + i) 



in 



det(^ M + m 2 ) 
{det(-df t) + mrf 



In 



+ In 



det (?{o,i +m 2 ) 
det(-<9 ( 2 =0) +to 2 

det(7£,_i l+ i 



m 



det(-d 



(<-§) 



+ m 1 



(2.38) 



+ ln 



det(^ + i |f _i +m 2 ) 
det(-<9 ( 2 /+1) + m 2 ) 



In 



det(76, t+ i + m 2 
det(— d 2 ^ + m 2 ) 



+ In 



l+ 2 ,l 2 



+ to 



det(H 
det(-d 2 +m 



while r^^A; to), the 'high' partial wave contribution, can be expressed by the form 



1>L 



(A;m) 



25tl{jTT(^-^)jr* ^ 
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~l J^f ln $ I d * x ( tr {FM + l tT {F ^ F ^) } ' (2 - 39) 



where we defined 



ft(r, r; s) = (21 + l) 2 {tr G M (r, r; s) + G?,_i >I+ i (r, r; s) + Gj+i^i (r, r; s) 
-2Gf ee (r,r; s) - G£|(r,r; s) - Gf_ e |(r,r; s)} - {G M+1 (r,r; s) 

+^ + i, z _i(r, r; S ) - Gf»(r, r; s) - G* e ! (r, r; *)}, (j = L + i, L + 1, • • •) . (2.40) 



In (2.38) and (2.39), various partial wave contributions have been grouped (with correct 
degeneracy factors) in such a way that we may have both the small-r and large-r behaviors 
of the effective potential matched when the radial differential operators figuring in within 
each group are taken together. (This grouping is especially important when we discuss 



the massless limit of the effective action (see Part A of Sec. IV)). In Case II background 
with (3 < (and finite R), the above partial- wave grouping is not really required; but it is 
desirable to have a procedure applicable to all cases we will consider. 



We remark that (2.37), even without taking the limit L — > oo, corresponds to an exact 
relation. But, only by taking L to be relatively large, this relation can be put to use as 
a powerful calculational tool for the effective action — this is because, for large L, two 
independent means to determine separately the high and low partial wave contributions 
become available. This is the key element of the partial-wave cutoff method. In subsequent 
sections, we shall obtain the expression for T^^A; m) as a ^-series form analytically while 
the other piece, T^j^A; m), in our chosen background fields will be found numerically (and 
semi-analytically in the small mass limit). If the sum of the two contributions yields an L- 
independent result, it is of course the sign that we have secured exact result for V^(A; m). 



III. HIGH PARTIAL-WAVE CONTRIBUTION TO THE EFFECTIVE ACTION 



In this section we shall calculate the large-L form of the expression (2.39) to desired 
accuracy. To that end we need a systematic large-/ approximation to the function Qi(r, r; s), 
which is valid uniformly for s in the range < si 2 < 0(1) pQ. For this one might consider 
using the WKB approximation [3J [15]; but, with a matrix- valued potential, it is not trivial 
to apply the WKB method directly. Alternative methods, yet serving our purpose in a 
satisfactory way, were found in [H [16]. These latter approaches are not only simpler but 
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also easily extendable to the case involving matrix-valued potentials. In this paper we will 
specifically use the method of [IB] to generate the desired y-expansion of Gij(r, r; s), because 
of its convenience in dealing with a matrix-valued potential. 

The idea of [16] is to rewrite the proper-time Green function, introducing the momentum- 
like variable p, in the form 



Gi,j (r, r'; s) 



dp I 

— (f 
2tt \ 



-sH 



1,3 



p) (p \ r ') 



& e ->[-%+Vu(r)] e ~ip(r-r') 



2tt 



(3.1) 



and to find a convenient way to study directly the r' — r limit of this function. Then, 



after moving the last Fourier factor e %vtyT r '- ) in (3.1) to the left of the differential operator 

r)r' 

G t j (r, r; s) 



d r (= J^), we may take the coincidence limit r' — r to obtain the following representation 



-oo 2tt 
00 dp 

-oo 2n 



K(r,p; s). 



(3.2) 



The function K(r,p; s) introduced here can be identified with the solution of the differential 
equation 



\ds 

under the boundary condition 



{d r -ipf + V Uj {r) \K(r,p; s) = 



K(r,p;s = 0) = I. 



(3.3) 



(3.4) 



Because of the connection (3.2), the desired large-/ series for Gij(r,r; s) will follow im- 



mediately if we have an appropriate development for the function K(r,p; s). For our inves- 
tigation, it is convenient to regard rescaled variables t = si 2 and q = f (rather than s and 
p) as independent variables of the function of K. Also we express the potential Vij(r) (see 



(2.36)) as a series in j, i.e., 



1.-. 



1 .-. 



V iA r ) = l ^ + 1 V ^ r ) + W + jVsir) + -Vt(r) + 



P 



(3.5) 



where Vfe(r) are some /-independent (matrix) functions. Then (3.3) can be written in the 
form 



d_ 
Of 



l9 --i 1+ — 
I dr I r 2 



VAr) , V 2 (r 



I 



+ 



P 



+ ■ ■ 



K[r, ql; p 



0. 



(3.6) 
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This equation may be solved by positing the ansatz 

K[r,ql;p \ =K (r,q,t) <1+ + - 2 + ■■■>, (3.7) 

where the function K Q (r, q, t), the I — > oo form of K(r, ql; 4), is chosen to satisfy 

(wt +q2 + ^) Ko{r > q > t) = - (3 - 8) 



The solution to (3.8), subject to the boundary condition (3.4), is 



K (r,q,t) = e-<« 2+ ^. 



(3.9) 



To determine cik(r,q,t) (k = 1,2, •••), we may plug in (3.7) and (3.9) into (3.6) to obtain 
following recurrence relations for them: 

2 

-n, (A; = 1,2,..-) (3.10) 



da k 



(d 8t\ (d 8tV * - 

-2iq — + — a fc _i + — + — a fc _ 2 - > V/„a fc - 



with a = 1 and a_i = 0. Then we can find a^'s successively (with a,k(r,q,t = 0) = for 
all k — 1, 2, • • •): here, the first few terms which are indispensable for our calculation using 



(2.37) are 



ai = —tV\ — 



8iqf 



a 2 



a 3 



iqt 3 - t 2 - 32g 2 t 4 16g^t 3 64t 3 12t 



r 



q 1 ^ m A 



32q 2 t 5 16<ft 4 64t 4 12t 



.2+4 



y« ■ J 



3r 6 



?' 



2 



(3.11) 
(3.12) 



2g 2 t 3 



3r 6 



+ 



-^r^i 2 + - ^ (2 W + t^vi) + - ^v? + t (v^ + w) 

tV + 256 ^ 3t6 128 ^ 5 + 32zg 3 t 4 512zgt 5 + 25Qiqt 4 64iqt 3 



3r 9 r 7 

512g 4 t 8 512g 4 t 7 2048g 2 t 7 , 384g 4 t 6 , 6016g 2 t 6 , 2048t 6 , 4zgt 5 T - 3 



3r i2 



.10 



3r 12 



3r 9 
+ 



(3.13) 



q 2 t 5 32t b 6t 4 lQq z f 



+ 



2 + 6 N 



+ 



8g 2 t 5 



+ 



+ 



x , 3r 6 r 4 

3328t 5 /t 3 g 2 t 4 ^ 
5r 10 + I 3"" IT 



4z'gt z 



9r 12 3r 3 

64g 4 t 5 1472g 2 t 5 

r J V 6 



(vm + V 2 V 1 ) + ~ + 2V{VA 



2,4 (ViVl' Vi'VA 4 , Wi \ , t 4 f/ 4 , 240g 2 t 4 488t< 

g M 2 + 6 +q 4 + 12 + 6 + 24^ 1+ r« + r* 



z'gt 3 
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~ (v?v 2 + v 2 v 2 + v 1 v 2 v l ) + l - (v 2 + v.Vs + vm) 

(64iqt 4 25Qiq 3 t 7 , 128iqH 6 , 5l2iqt 6 32iqH 5 256iqt 5 ^ 
I r 5 



+ 



+ 



3r 9 r 7 3r 9 

'32g 2 t 5 _ 16g 2 t 4 _ 64t 4 12t 3 \ - 

^6 ^4 3 r 6 ^ r 4~ J 2 ~'~ 

'8g 2 t 4 16t 3 \ [32iqt 4 16iqH 5 



'2q 2 t 3 t 2S 
3 2 / 

) V 1 + iqt 2 V 3 ' 



V" 

V 2 



'I6ig 3 t 5 32iqH 6 64iqt 5 32iqt v 



3r 3 



^ 3r 3 



3r 3 



^"6 /y 

1.3 „-„3+4 



6 



2igt 3 zg 3 r 



(3) 



(3.14) 



where VI ,V fc " and denote the first, second and n-th derivatives of Vk(r), respectively. 



We may use the expansion (3.7) for the function K in (3.2) and carry out the p-integration. 
Since 



dp 
2^ 



K (r,q 



V 



t = sl' 



TVS 



(3.15) 



the result is the following j-expansion of the quantity Gy(r, r; s) (which can be used even 
when si 2 ~ 0(1)): 



-sl' 2 - 



Gi,j(r,r; s 
1 



+ 



I 3 



y/Ans 
4(s/ 2 ^ 3 



1 + 



1 r 



-sZ 2 Vi 



+ 



P 



3r 6 



2\4 N 



3r 3 



4(s/ 2 ) 3 16(s/ 2 ) 



3r 6 



(sl 



2 ) 3 -o (Sl 2 ) 2 f- - - - 

+ [VM + V 2 V 1 



6 



V? - sl 2 V 3 



+ 



I* 



2(sl 



3r 3 



5- (v.v; + v[% 



\(sl 



2\5 



2(sl 



2\4* 



3r 6 



A{sl 2 f vij {sl 2 f 



3r 3 



6 



VxV{' Vl'Vi -o- - -o 

+ + - V 2 V 2 - V 2 V 2 - V 1 V 2 V 1 



(sl 



2\4 



+ 



/4(s/ 2 ) 3 16(sZ 



2\4 K 



3r 6 



^2 + 



(sZ 



2\2 



K 2 + W 3 + V 3 Vi - \ 



- sl 2 V 4 



24 
128(sZ 2 ) 6 



9r 12 



704(sZ 2 ) 5 40(sZ 2 ) 4 



5(sZ 



2\3 



15r io 



+ 



(3.16) 



The (matrix) functions Vi(r), V^(r), ■ ■ ■ here should be found through (3.5), for the potential 



Vy given by (2.26)-(2.29); hence, Vi(r),V 2 (r), ■ ■ ■ can be expressed in terms of our profile 



function f(r) (see (1.1)). 



Using the form (3.16) with (2.40), the systematic large-/ approximation of the quantity 



Qi(r, r; s) can also be obtained. We may then use the result together with (2.39) to evaluate 
r^(A; m), according to the general procedure we detailed already in [IJ|2j[T6]- That is, the 
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/-sum appearing can be performed with the help of the Euler-Maclaurin summation formula 
and this is followed by the integration over the proper-time variable s. One can then verify 
that the A — » oo limit is indeed well-defined, for the potentially-divergent proper-time inte- 
gral gets canceled by the renormalization counterterm contribution. After these, somewhat 
lengthy but straightforward, calculations, we have found that the quantity T^^A; m) in 



the gauge background (1.1) can be expressed as a ^-series of the form 

'Q 2 (r)L 2 + Q 1 (r)L 



1>L 



(A;m) 



dr 



+Q log (r) ~ 



+ Qo(r) + 



(3.17) 



where u = dl 
4r 



2L 



Qo(r) 



r 
" 2~ 
r 

24m 7 



Q 2 (r) = -f(r) (r 2 f(r) - 
u v 

Q 1 (r) = ^ (Su 2 - l) /(r) (r 2 /(r) - l) , 

Qiog(r) = V {/(r) 2 (6r 3 /(r) - 20) - r 2 f{r f - 10r/(r)/'(r) - 4r 4 /(r) 4 + 20r 2 /(r) 3 } , 

r 3r 2 /(r) 2 (-24rV/'(r) + 128m 6 - 55m 4 + 30w 2 + 5) (3.18) 
+2ru 2 {-ru 2 f"(r) + r 3 u 2 (Qu 2 + l) f(r) 2 - 3 (2m 4 + u 2 + l) f(r)} 
+/(r) {4r V/"(r) + 4r 3 (30w 4 + 5m 2 + 3) u 2 f'{r) - 3 (56?i 6 - 57u 4 + 24u 2 + 5) } 



+2rV (24u 4 - 13m 2 + 3) f{r f - 4rV (60u 4 - 7m 2 + 3) f(r) 

1 (r/ P ) 2 a 



If the formula (3.17) is used to our Case I where f(r) 



after the r-integration, the following result: 

2L 2 3L 

a a 



In 



a 1 
6 + 6a + 2 



r 2 (r/p) 2 «+l 

1\ 5a 



36 



(a > 1), one obtains, 



^ +0 Q) <3.19) 



For the mathematical validity of our formula (2.37), we do not need to know, in the right 
hand side of (3.17), the explicit form of the 0( j-) term or beyond, i.e., Q^i(r)jr-\-Q_ 2 (r)-^ + 



• - -. But, as noted in [TJ 121 HE], these -^-suppressed terms can be important to accelerate 
the convergence of our calculational scheme if the low partial wave contribution has to be 
evaluated by numerical methods. They can of course be found by keeping further higher- 



order terms in (3.7) and (3.16) above. The end results turned out to be quite lengthy: see 



Appendix |A| for the explicit expressions of Q_i(r) and Q^{t)- (These results are utilized 
in Sec. IV). 
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IV. LOW PARTIAL- WAVE CONTRIBUTION AND THE FULL EFFECTIVE AC- 
TION 



Our next task is to evaluate the low partial-wave contribution r^(A; m), given in (2.38) 



for the background field of the form (2.30) with certain specific f(r). By combining this 
with the high partial- wave contribution (calculated already in Sec. Ill) a la ( |2.37[ ), we can 



determine T^(A;m); then, by (2.31), the full fermion effective action T ren (A;m) follows at 



once. We use the Gel'fand-Yaglom (GY) method [12] to determine the ratio of two functional 



determinants in (2.38). Applying this method to functional determinants involving ordinary, 



non-matrix-type, differential operators is well-known, and we can thus write [2] 



det(Hu + m 2 



lim 



1pl,j(Re 
'.free/' 



(ifj^l, or/ = j = 0), 



det(-<9 { 2 + m 2 ) R e ^oo ?/f ee (i? e ) ' 
where 4>i,j{ r ) an d ipf ee (r) denote the solutions to the radial differential equations 



(4.1) 



{Ui,j + m 2 )^i,j{r) = 0, 



(4.2) 



c) 2 



m 2 )ip t 



free / 







with following small-r limit behaviors 



r — > : ipij(r) ~ 1 • r 
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^ free (r) ~ 1 ■ r 



2/ 



(4.3) 



But, in (2.38), there are also functional determinants involving 2x2 matrix differential 
operators, i.e., Hu + m 2 . For these functional determinants, we have to use the generalized 



GY formula [Tf] 



det{H ltl + m 2 



lim 



det(ipi jCe p(R e )) 



{det(-d 2 )+ m 2 )y ^ ^(Ref 
where ipi^pij") (a, {3 = 1,2) denote the solutions to the differential equations 

CHi,i + m 2 ) a ^i^x(r) = 0, 
(Hij, + m 2 ) Q/3 ^,/32(r) = 

with following small-r limit behaviors 



(4.4) 



(4.5) 



r 



ipi,2i( r ) ^,22 (r) 




(4.6) 
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Note that the two equations in (4.5) may be written as a single matrix differential equation 



{n^i + m 2 ) ¥,(r) = 0, 



^,11 (0 ^,12(0 
^,21 (r) ^,22 (V) 



(4.7) 



With nonzero mass m and a generic radial function f(r), exact forms of the GY wave 
functions ipij(r) and *&i(r) are usually not available. But, in the massless limit, a certain 
analytic procedure can be developed to find these wave functions. In fact, with the back- 
ground field of our Case I, complete GY wave functions may be obtained by this procedure. 
So, below, we shall show how our general procedure can be applied to calculate the effective 
action in our chosen backgrounds in the massless limit first (Part A). This will then be 
followed by the corresponding discussion with m ^ (Part B), which requires extensive 
numerical analysis. 

We note that, in performing this effective action calculation, it is convenient to set the 
length parameter p (entering our background fields), and sometimes the normalization scale 
p also, to be equal to 1. This does not amount to a loss of generality. It is related to the 
fact that, from simple dimensional argument and the way the normalization scale p enters 



r re n(^4; fra) (see (2.4)), the modified effective action T(A;m) in our backgrounds, defined by 
the relation 

T mn (A; m) = 2 - ln(//p) j -0^ tr(F^F^) + T(A; m), (4.8) 
should be a function of dimensionless parameters not involving p, i.e., a function of mp and 



a for Case I, and a function of mp, R/p and j3 for Case II. According to (4.8), we now have 
that 

T rcn (A; m)\ p=fM=1 = f(A; m)\ p=1 , (4.9) 

viz. by calculating the effective action with p = p = 1, we have calculated T(A;m)\ p= i. 
But the modified effective action T(A;m) for arbitrary p-value follows from T(A;m)\ p=1 by 
dimensional considerations — just regard the numbers assumed for m and R in T(A; m)\ p= i 
as denoting the values of mp and R/p, respectively. With the quantity T(A; m) thus found, 
the corresponding effective action r ren (^4; m) for arbitrary values of p and p is provided 



through (4.8) 
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A. Fermion effective action in the massless limit 



For m = the above GY equations exhibit a special feature of factorizability. To show 
this, we note that the differential operator representing — D 2 — |^l^c a -^v in the radially 



symmetric background (2.30) (see (2.16)) can in fact be decomposed as the product of two 
linear differential operators, according to following, directly verifiable, relation 



2 



3 d 



+ — L 2 + 3r 2 f(r) 2 + 8/(r)f • + 4 Uf(r) + r/'(r) - 2r 2 f{ 



S-T 



= -(■?- + - + • S + Arf(r)S- f}(^-- • S - Arf(r)S ■ f) . (4.10) 

\or r r J \or r J 

Our operator Hij is nothing but the restriction of this operator to the partial waves with 
quantum numbers I and j. For given values of I and j, S ■ T is represented by the form in 



(2.25); similarly, for ■ S (= \{J 2 — Q 2 — S 2 — 2S ■ T)), we have the representation 
• S - 



I 
2 



l + l 



( 



2(21 + 1) 



1(21 + 3) 2yjl(l + 1) 



^ 2- x /Z(Z + 1) (Z + l)(-2/ + l) 







if j = 1 + 1, 

if j = 1-1, 

if j = 1^0, 

if j = 1 = 0. 



(4.11) 



Also a remark as regards the j = I = partial wave: for our backgrounds having a nonzero 



Pontryagin index, there will be a normalizable zero mode of the operator (4.10) 



Based on the above observation, we can recast the massless GY equation Hijip(r) = 



appropriate to j = I + 1 (and hence S ■ T — > ~ and Z/ + ) ■ S — > ~) 



as 



dr r r j \dr r 



rf(r) ip(r) = 0. 



(4.12) 



Therefore, for the GY wave function, we may well look for the solution to the first-order 
equation 



' d 21 
dr r 



rf(r) ip(r) = 0. 



(4.13) 



This way, the GY wave function with the correct small-r behavior is obtained: 

i>u=i + i(r)=r 2l J» r ^ dr \ 



(4.14) 
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As for the GY equation with j — I — 1 (and so S ■ T — > |, ■ S — > — ^jr), i.e., for the 
equation 

d 2(/ + l) 



9 3 2(/ + l) 

7T + - + r/ r 



+ 



- rf(r) ij}[r) = 0, 



Qr r 

the situation is not quite the same. In the latter case, solving the first order equation 

(d 2(1 + 1) 



(4.15) 



dr 



rf(r) ipi(r) = 



results in a solution of the form 



^( r ) =r -2a+l) e / o r n/(r 1 )dr 1 



(4.16) 



(4.17) 



which is singular as r — » 0. For GY wave function we thus have to look for another 



kind of solution to the second-order equation (4.15). For such solution ip2(r), we here put 



^(t) = ip\(r)a(r) and then use the equation (4.15) (together with (4.16) for ^i(r)) to obtain 

(d 3 2(1 + 1) 



dr^r 



_l_ rf(r) ifji(r)a'(r) = 0. 



By solving this equation we are led to the expression 

a (r) = r e J o LJ L . 



(4.18) 



(4.19) 



and this in turn lead to the following form for the second solution fair) = ipi(r)a(r): 

^ 2 (r) = r-^+Vo^ 1 ^ 1 ^ 1 f rf +1 e- 2 ^ 2r ^ dr ^dr 2 . (4.20) 



21 

For small r, this function behaves as fair) ~ 2(21+1) • Hence we can identify the appropriate 
GY wave function with 2(2/ + 1) times this function, i.e., 

^, i= ,_i(r) = 2(2/ + i) r -2P+i) e /o P ^/(^)^ f r^e' 2 ^ Tlf ^ )dr Hr 2 . (4.21) 



Massless functional determinants for partial waves corresponding to j = I ± 1 can be 
evaluated using the GY wave functions in (4.14) and ( 4.21[ ). Especially, for our Case I, i.e., 



fi r ) — \H(r) with the function H(r) as given in (2.21) (while taking p = 1), the exact 



wave functions are given in terms of hypergeometric functions: 
^ J=l+1 (r) = r 21 (r 2 « + l)^, 



1 2/ + 1 2/ + 1 



Aj=i-i(r) = r 21 (r 2a + 1 2a 2 F X [ -, — ; — + 1; -r 



2a 



a a 



a 



(4.22) 
(4.23) 
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Since we have ip\ 



free / 



r 21 with m = 0, we may now readily calculate the asymptotic wave 



function ratios in ( |4.1 ) to conclude that 
In 

In 



det T-Lij = i + i 


~ In 


ipi d= i+i(R e ) 


_ det(-ag)) _ 


[ ^r{Re) \ 


det Hi,j=i-i 


~ In 


ll>l,j=l-l{Re) 


_ det(-«9 ( 2 ) _ 





In R P 



In R e + In 



2/ 



(4.24) 
(4.25) 



viz., the corresponding functional determinants individually are not well-defined. This prob- 
lem, noticed in a similar context also in Ref. j3], occurred because of our setting m to be 
exactly zero. As a matter of fact, the asymptotic ratios found with m set to zero are in 
general not the same as the massless limits of the asymptotic ratios calculated assuming 
nonzero mass. (We elaborate on this aspect in Appendix [Bj . But, if one makes a 'good' 
grouping of different partial- wave contributions, the two results for the group coincide [3]. 
In our case, such good grouping is provided by the way we combined various partial-wave 



contributions in (2.38). (This is justified in Appendix pi). With this understanding we may 



apply our results (4.24) and (4.25) to the particular combinations entering (2.38), to write 

det Hj.i 



hi 



det H x _ 



det 



-f) 2 



+ In 



det 



-f) 2 



In! 



21 + 2 
.27+1. 



I = 

2' ' 



(4.26) 



and 



In 



det Hi 



,3=1+1 



det(-<9 ( 2 ) 



+ In 



det ft 



i+y=i- 



det 



-B 2 



In 



'2/ + 2 x 
21 + 1 



(4.27) 



[We here remark that, although (4.27) follows also from considering the small mass limit, 



(4.26) does not; but, if the formula (4.26) is used together with (4.50) below, the value ob- 
tained for the total sum becomes also consistent with the massless limit of the corresponding 
massive expression]. Analogous considerations may be given to our Case II as well. But, to 



obtain the corresponding values for the quantity in the right hand side of (4.26) or (4.27), 
numerical integration will be required. 

Our next task is to study the functional determinant from the j = I = partial wave. In 
this case, S ■ T — > — | and • S — > and so we have the GY equation 

( d 3 _ „, A / d 



dr 



3r/(r) 



dr 



+ 3r/(r) ^(0 = 0- 



(4.28) 
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The GY wave function may then be identified with the solution to the first order equation, 
i.e., 



j=l=0 



-3 J Q rif(n)dn 



(4.29) 



But, if the function H(r) = r 2 f(r) approaches 1 as r — > oo (i.e., for Case I and also for Case 
II with (3 > 0), we find that V ; i( r )|j=«=o ~ as r — ► °°- This corresponds to a normalizable 
zero mode mentioned earlier, and the related functional determinant vanishes. In this case 
our interest will naturally be in the expression when a small mass m is included; in the GY 
approach, this requires the knowledge on the asymptotic behavior of the GY wave function 
satisfying the equation ("Ho,o + m 2 )^ = with small but nonzero m. The latter GY wave 
function, which we denote as V ; i=«=o( r ), can also be constructed using the method of Ref. 



[I]. Based on such analysis (see Appendix |B| especially (B13)-(B20)), we then obtain, say, 
for our Case I, the following result 



In 



det (7^o,o + m 2 ) 
det(-dg =0) +m 



In 



l/jj=l=o(Re 



In m + In 



r i 



2r © 



(4.30) 



when m is small. As this corresponds to the very first term of the first group in (2.38) 



we may combine this result with the small-mass-limit form of the second term in the same 



group. The relevant result, for Case I, is (this follows from (Bll)) 



In 



^=0,j=l (-Re) 



Hence, in the small-m limit, we have 

det(U 0fi + m 2 ) 



In 



+ m z 



det(-df l=0) 
= 2 In m + In 



In 



lnm — In 4. 



det (?{o,i +m 2 ) 
det(-d 2 +m 2 



(4.31) 



8T(S) 



(Case I). 



(4.32) 



[For our Case II with /3 < 0, the GY wave function (4.29) has a nonzero asymptotic value 
(i.e., lim^oo ipi(R) ^ 0) and this limit value determines the j — I — massless functional 
determinant] . 

Let us now turn to the case j — I ^ 0, i.e., the case where 2x2 matrix differential 



equation (4.7) is relevant. Here, in dealing with the boundary condition (4.6), the non 



diagonal matrix form given for L(+) ■ S in kill is not very convenient. Therefore, we 
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perform a unitary transformation, ■ S — > U(L^ ■ S)W , with 

(UW = I 



U 



2Jl(l + l) 



21 + 1 ^ 2 V /Z(Z + 1) 1 
to find the following diagonal form for ■ S: 



(4.33) 



L<+> • 5 



1(1 
— Z — 1 



(4.34) 



Under this unitary transformation, 5* ■ T (originally given by the form in (2.25)) remains 
unchanged, i.e., 

(4.35) 



S ■ T 



4(2/ + 1) 



21 - 3 Ayjl(l + I) 
v 4 V //(Z + 1) -2/ + 1 



Then, based on the factorized form (4.10), we may first consider the first-order matrix 
equation 

/ 



d 2 11 
dr r I o -I - 1 



2/ + ] 



-2/ - 3 4 J/(/ + 1) 



^ 4 V / Z(Z + 1) -2Z + 1 



*(r) 



0. 



(4.36) 



where ^(r) is a 2 x 2 matrix (see (4.7)). For a general radial function f(r) (assumed to 



be finite for r — > 0), it will be unwieldy to exhibit the solution to this matrix equation in 



an explicit manner. So just let a 2 x 2 matrix function ^i(r) denote the solution to (4.36) 
which has the following small-r behavior 

' 1 • r 21 ■ r~ 2l ~ 2 



r -> : VfMr) ~ 



0-r 



21 



1 • r 



-2Z-2 



(4.37) 



A comment as regards (4.37) might be appropriate here. Note that, if (4.36) were regarded 



as an equation for a 'column vector', our '2x2 matrix' ^ would comprise two column 



vector solutions to (4.36). Then (4.37) is equivalent to the statement that we require two 
independent solutions to the column vector equation (4.36), say ^^(r) and \l/^(r), having 



the small-r behaviors 



r : ^S 1} (r) ~ 



.21 



] r -2/-2 



(4.38) 



Notice that the solution to the first-order equation (4.36) is not the one satisfying our 



boundary condition (4.6) (with the above unitary transformation taken into account). Then 
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the full second-order GY equation for j = I ^ should admit a different kind of solutions. 
To find such solution SE^r), we put ^O") — ^ , i(r)A(r) (A(r) is a 2 x 2 matrix function to 
be determined) and use the form with the second-order equation 

[ ^ + 3 + 4 L (+) ■ S + 4rf(r)S- T] I ~ - 4 I (+) • S - 4rf(r)S ■ f) *x(r) A(r) = 0(4.39) 
\or r r ) \or r J 



for • S and S ■ T given by the matrices in (4.34) and (4.35). This then reduces to the 
first-order equation for <fr(r) = ^ r 1 (r)A'(r): 



9 3 2 
d~r + r + r 



( 



I 
-/ - 1 



+ 



rf(r) 
21 + 1 



( 



-21 - 3 4 J 1(1 + 1 



V 



4 v //(/ + l) -21 + 1 



$(r) = 0. (4.40) 



If «&3(r) denotes the solution to this equation with the small-r behavior 



r — > : $3(r) 



■y . y,— 2Z— 3 g . j,2i— 1 
Q . 2Z— 3 j . j.21—1 



(4.41) 



the desired second solution ^(j") can be identified with 

/r 
[*i(ri)] _1 * 3 (ri)dri. 



(4.42) 



Thanks to (4.37) and (4.41), this second solution has the small-r behavior 



: Voir) ~ 



■1 • r 



-2Z-2 



0-r 



21 



2(21 + 1) I Q . r -2Z~2 ^ . r 2/ 



(4.43) 



Now, for the solution satisfying the GY boundary condition (4.6), we consider a linear 



superposition, i.e., &i(r) = ^i(r)Ci + ^/ 2 ( r )C2 where Ci and C 2 are suitable constant 
2x2 matrices. By this consideration, we can identify the GY wave function for j = I ^ 
with the expression 



*i(r) = #i(r) 



1 




2(2Z + l)* 2 (r) 




1 



(4.44) 



(in the basis where L^' ■ S and S • T are given as in (4.34) and (4.35)). Note that, if two 
column vectors ty[ (r) and ty^\r) (^^(r) an d ^^i 7 ")) are used to represent our solution 



^i(r) (^(r)), (4.44) identifies the desired GY solution ^/(r) with the matrix formed by 
two column vectors \&i (r) and 2(2/ + l)^^\r). 
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For our Case I the above matrix functions ^i(r) and $3(7") can be found explicitly: if 
\E r i(r) nm ($ 3 (r) nm ) denotes the nth column and mth row of ^ , 1 (r)($ 3 (r)), we have 



*i(r)n 
*i(r)ia 
*i(r) 21 



22 



3v;i2 



$s(r 



3VJ21 



31' j22 



r 2* / 2. + ^"i f _ 1 » 2J + 1 _ r2Q 
v y V a a a 



2^1(1 + 1) (r 2a + 1) 2 - r 
(2/ + !)(-« + 2/ + 1) 



„2a-2«-2 



ft 



ft 



ft' 



Zr 2Z + 



2o 



'2/ 



ft ft ft 



ft ft ft j J 

r -a(i+i) ( r 2a + ^-^ r 2 ( a _ !)( 2 / + l)(- a + 2/ + 2) (r 2a + 1) r 2a 



(4.45) 



2(2/ + l)(-ft + 2/ + l) [ 

/ 1 21 + 2 2Z + 1 

x 2 Fi 2--, 2 ;3 ;-/■ 

V ft ft ft 



-2ft + 2/ + 1 

2 " ^ + (2(2/ + l)(-ft + 2/ + l) 



-2[ft + 2l(a - 21 - 3) - l}r 2a 1 2 F a ( 1 - -, 1 



21 + 2 



21 + 1 



,2a 



-r 



ft 



and 

*3M„ = r" 2 '- 3 



2^/(Z + l)(r 2a + l)^r 
(2/ + l)(ft + 2/ + 1) 
/r -2i-3 ( r 2« + ^ 



ft ft 



2^ „2a+2Z-l 



ft 



/ 1 21 + 2 21 + 1 9 „ N 
2F1 1 + -, 1 + 2 + -r 2a 



ft: 



ft 



a 



I 2a 



sjKi + 1) 



r 2Q + l) 2 Fxil + -,l 



21 



;1 



2/ + 1 



„2a 



ft 



ft 



-2 -Pi I 



2/ 21 + 1 



5 ? 
ft ft 



-; -r 



2a 



ft 



(4.46) 



,21—1 ( m 2a 



r 2a + j_J 



2 a 



2(2/ + l)(ft + 2/ + l) 



(2(2/ + l)(ft + 21 + 1) + 2 [a + 2Z(a + 2/ + 3) + l]i 



.2a 



2/ + 2 2/ + 1 



x 2 F 1 \l + -,l + ;2 + 

V ft ft 



,2a 



ft 



2(ft + 1)(2/ + l)(ft + 21 + 2)r 2 « (r 2 « + 1) L + i 2 + 21 + 2 g + 21 + 1 . _ r2a \ 1 
2ft + 2Z + l y a' a' ft' yj' 

The second matrix solution ^ 2 { r )j and also the GY wave function Vf^(r), will be given 



using these results. But, for the determinant ratio (4.4), all we need to know is the form of 



det *f?i(r) for large r = R e , and, because of (4.42) and (4.44), it can be recast as 



det¥,(fle) = 2(2/ + 1) (det*i(i2e)) ([^linT^siri))^ dr x . (4.47) 
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Here, det *&i(R) is easily calculated (without using the explicit forms in (4.45)) if one notes 



that, for ^i(r) satisfying the first order equation (4.36), 



ln[det(r*i(r))] = trln(r*i(r)) 



tr 



tr 



+ *;(r)*i(r) 



12/ 



-i 



dr 



| rf(r) "2/ -3 4^ + 1 
-I - 1 j 21 + 1 1 + 1) -21 + 1 



rr 

—2 / rf(r)dr. 
Jo 



Hence, for f(r) 



and accordingly 



r 2 (l+r 2a ) ■ 



ln[det(r*i(r))l = -- lnfl + 



,2a 



det*a(i^ e 



1 • 



On the other hand, from (4.45) and (4.46), we find 



r -> oo : (*i(r) x $3(r^ 



22 



(2/ + i)V +3 r(^) 4 
4/(/ + i)rg) 2 r(^) : 



Using these results in (4.47) gives rise to 



(21 + lfRfT 



3 D4«r / 2Z+1 



8/(/ + i)^r(f) r(^)' 



Based on this, we have 

det "Hu 



In 



{det(-%)} ; 



. det*,(i2e) . 

In — t — - — — = in 

^(^e) 2 



(2/ + i) 3 r(^) 4 
8/(/ + i)2r(f) 2 r(^) : 



dr 



(4.48) 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



In (2.38), this result may be used in conjunction with that in (4.26). 

To obtain the quantity T^ L (A;m) in the small-mass limit, all that is needed now is to 
consider the sum of various partial-wave functional determinants discussed above. For the 



backgrounds corresponding to our Case I, we find from (4.32), (4.50), (4.26) and (4.27) the 
following result for the sum: 



T^jj (A] m) = — In m — — In 



1 



~\ E ^(2/ + l) 2 ln 



rfi + ijr 



a I \ a 



(2/ + i) 2 r( 2 ^)^ 



4/(/ + i)r(f) r 
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2/+2 



+ ln 



'21 + I s 
21 + 2 



(4.54) 



For large enough L this quantity can be computed as follows. Here notice that the quantity 



inside the curly brackets, waiting for the Z-sum in (4.54), can be approximated for large I by 



4| _ 2 

a a 



a 2 +3 
6la 



+ Oi^j^j. Therefore, if L is large enough, we obtain from (4.54) 

1 



r 



(-) 

KL 



2L 2 3L n , r . , fa 1 
a a V6 6a 2 



— In m In 

2 



r(l+M|) 
8r (t) 



+ C(a)+0[-), 



(4.55) 



C(a) being given by 



C(a) = — £ ^(2/ + l) 2 ln 



(2i + i) 2 r 



2«+l 



41(1 



+ ln| 



'2Z + 1 
21 + 2 



4/ 2 a 2 + 3a + 1 

+ - + - + cl 

a a ota 



(4.56) 



The constant C(a) may be evaluated numerically. Observe that 0(L 2 ), O(L) and O(lnL) 



terms in (4.55) match precisely those of the high partial-wave contribution given in (3.19) 



We thus obtain the unambiguous result for their sum, i.e., for the quantity T^ e 2(A;m) 



according to (2.37): 



ren l^J 



m) 



-ih. 

2 



r 




r 


id 


sr 







a 1 1\ 

— I 1 — -mm 

6 6a 2/ 



ha _ 47 
f 36 ~ 72a 



C(a) 



(4.57) 



Using the result (4.57) with (2.31 ) then provides us with the exact expression for the small- 



mass-limit form of the renormalized fermion effective action: i.e., for our Case I backgrounds, 



r ren = 2r(- ) + inQ 



a 2 + l 

- ln(mp) + 3 ^ ln(/xp) + C(a), 



(4.58) 



where we reinstated the p-dependences, and C(a) is given by 



C(«) = (7-ln2)|'f + j£ 



111 



r 




r 




sr 







+ 



5a 



47 



18 36a 



2C(a). (4.59) 



Note that the numerical factor 



multiplying ln(mp) in (4.58) reflects the value of the 



classical Yang-Mills action, i.e., / -^^^(F^F^) 



a 2 +l 
2a 



for our backgrounds. For a 



i.e., when the background field corresponds to a single instanton solution, we have C(l) 
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4£'(— 1) + + J; = —0.291747; then, from our expression (4.58), the previous calculation 
of 't Hooft [6] is recovered. Values of C(a) for different choices of a are given in TABLE 
[ij Especially, for large enough a (so that H(r) may resemble a step function 9{r — p)), one 



a 


1 


2 


3 


4 


5 


10 


20 


C(a) 


-0.291747 


-0.269189 


-0.378112 


-0.590437 


-0.883495 


-3.16105 


-10.0277 



TABLE I. Values of C(a) for various a. 



can derive from (4.56) an appropriate asymptotic formula for C(a); based on this, one gets 
the expression 



ln(mp) + 



a 



3a 



ln(pp) 



a In a fl 



ln2l . / 4\ ^/W 



(4.60) 



which proves to be quite accurate, say, if a > 8. Note that, for very large a (i.e., when H{r) 
is almost step-like), this one- loop result can be very large (the — a ^ term can exceed the 
classical action value); but then, higher-loop terms can also be significant. 

In our Case II backgrounds the functional form of Hu(r; R, f3) is such that no simple 
expression can be obtained (even with the simplification introduced above) for the small- 
mass-limit form of the effective action. Hence we shall be content here with exhibiting 
certain feature concerning the massless fermion effective action in our Case II backgrounds 
with (3 < 0, i.e., for the case of instanton-antiinstanton-type configurations shown in Fig. 
[ijd). Actually, for the present discussion, we may take the function H{r) = r 2 f(r), entering 



(4.61) 



the background field (2.30), to have the general form 

H(r)=H a (r){l + H b (r-R)}, 

with the functions H a (r) and Hb(r — R) broadly observing the patterns shown in Fig. [3j 
Then it will be of interest to know the behavior of the fermion one-loop effective action as 
R, the instanton-antiinstanton separation, becomes large. With finite fermion mass, one 
expects that it should be approximately equal to the sum of the individual contributions 
from the instanton and the antiinstanton. (This is also borne out in our numerical study, 
presented in Fart B). But, with negligible fermion mass, this is known to be generally not 
true [3 18] — there exists long-range interaction between the instanton and the antiinstanton. 
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H a (r) 
1 



Hdr-R) 



R 



R 



(a) (b) 
FIG. 3. Schematic forms of the functions H a (r) and Hb(r — R). 



We would like to identify such long-distance interaction term in the massless effective action 
when the background field has the above form. 

Since the Pontryagin index is zero for the above background, we have Y rcn (A\m) = 
2T^(A;m). Then note that, among various partial wave contributions to T^~2(A;m = 0), 
the I = j = partial wave term is rather special. If the background field have had only 



the instanton part (i.e., without the if^ part is (4.61)), a normalizable zero mode would 
have been present in this partial wave, and hence a divergent contribution to the effective 
action. But with the Hb part included (i.e., if an antiinstanton is also present at some 
distance r = R), there is no normalizable zero mode in any partial wave term. For the 
I = j = partial wave contribution, this amounts to a big change, from a divergent result 
to a finite one. If the instanton-antiinstanton separation R becomes quite large, we must 
then be able to see some, strongly independent, term (representing instanton-antiinstanton 
interaction at large distance) from this partial wave contribution. Further, when the mass 
value is sufficiently small, our numerical study (presented in Part B of this section) shows 
very clearly that the contributions from other parts do not generate significant long-range 
interaction term. Therefore, to extract the very long-distance interaction term, we may 
concentrate our study to a specific group containing the I = j = partial wave contribution 



of T^(A;m), i.e., according to our grouping made in (2.38), to that consisting of the 
I = j = and (I = 0, j = 1) partial wave contributions. 

With m = the I = j = and (I — 0,j — 1) partial wave GY wave functions, ?/>o,o( r ) 
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and ipo t i(r), are given by (4.29) and (4.14): 



3 Jo 



r H( ri ) 



(4.62) 
(4.63) 



Therefore the lowest angular momentum part of our effective action expression (2.38) be- 
comes 




det "Hq, 



o 



det(-<9 ( 2 



In 



— lim In 

2 T— >00 



(1=0)) 

^o,o (r) 



det H 



o.i 



det(-d, 2 



(1=0) J 



In 



det % 



o.o 



det Hq i 



^0,1 ( r ) 



r 00 H(r) , 
jo r 



(4.64) 



Now, with the form (4.61) for H(r), we may rewrite this quantity as 



H(r) 



dr 



l H a (r 



-dr + 2 



CO 



H a {r) 



+2 



or Ji r 

- H a (r)H b (r - R) + 1 



j o r 

/■A 1 

dr + 2 -dr. 
Ji r 



(4.65) 



We will take R to be large. Then the first two terms in the right hand side of (4.65 ) are finite 



and -R-independent. The third and fourth terms are independent but remains finite for large 
R. But we have also the last term, 2 ^dr = 21ni?, i.e., a term growing logarithmically 
with R. Based on this, we can now conclude that the massless effective action T ren (A; m) 
in the above well-separated instanton-antiinstanton background should contain a long-range 
logarithmic interaction term (of attractive nature), i.e., 

r ren ~ 4 In # + 0(1). (4.66) 

This is consistent with the observation of Refs. (7J |8] . 



B. Fermion effective action with m ^ 

With m^O, numerical integration should be considered to solve the GY equations. But, 
for a relatively large value of mass m, we have a totally different approximation scheme for 
the total effective action in the form of the large mass expansion. To acquire a measure on 
the validity of the latter scheme, we shall below summarize the appropriate formula of the 
large mass expansion first. Note that we here assume p — p, — 1; so without this assumption, 
m and R below become rap and — , respectively. 
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1. Large mass expansion 



One can obtain the large mass expansion of the one-loop effective action with the help 
of the Schwinger-DeWitt proper time expansion. Since the related procedure is described 
in detail in (TJ [2] for the case of scalar effective action, we will present only the final results 
that apply to our discussion with fermions. The large mass expansion for r^(A; m), in the 



background of the form (2.30), can be written as 



rH(A;m) 



In m + A {2) — + A {4) — + 
/i LM m m -t- /i LM 2 + ^lm a + 



(4.67) 



All of the coefficients in (4.67) involve radial integrals of certain polynomials of the function 



f(r) in (2.30) and its derivatives. Explicitly, they are of the following forms: 



A 



(o) 

LM 



a(2) 
^LM 



.4 



(4) 
LM 



2 Jo 

oo 

dr 



dr r 3 (/ 2 [20 - 6r 3 f ] + r 2 f' 2 + lOrff + Ar A f - 20r 2 / 3 







1 

40' 



r 4y(3) \-3rf + iQr 2 / 2 - 16/ 



1 

120' 



-80r 5 / 3 /" 



-3r 3 /" 2 + 510r 3 / 2 /" - 432r//" - 600r 6 / 4 /' + 60r 5 / 2 f 2 
+30r 4 /' 3 + H20r 4 / 3 /' + 540r 3 //' 2 + 930r 2 / 2 /' - 302rf 2 
-720//' + 120r 4 //7" - 139r 2 /7" + 400r 7 / 6 
-2160r 5 / 5 + 3120r 3 / 4 - 960r/ 3 ] } , 

r dr J _ /2 /(4) _ 113 4 /2 /(3) + li r 3 //( 4) + 99 2 (3) 
Jo I 5 ^ ^ 40 ^ ^ 20 ^ ^ 56 ^ ^ 

+ -r 7 f 2 f" 2 - 2r 5 ff" 2 + — r 3 / 2 /" + -r 8 f/' 3 - -r 7 /' 4 
15 JJ 40 ^ ^ 10 J J 30 J 

_ 124 r 6 ; ; /3 _ 61r 5 /2 ; /2 _ 173 r 4 f ,3 _ 2 r 7 /6 ^3 f _ 354 

15 60 5 L 



1GS0 . L 13r 2 (/ (3) ) 2 + 1675(f) 2 + 16r 2 /(V + 488r/( 3 )f 

157 r 4 /( 4) fl 299 3 (3) _ 4 (5) r_ 2 ?r2/2 _ n ; 
840 J J 210 J J 280 J L J J 

--r 6 f'f" 2 + -r 2 f'f" - -r 3 /f 2 [36r 4 f - 43 
40 J J 12 J J 60 J J L ^ 



r 5 / 5 25r 4 /" - 207r 3 f + 450 



15 
1 



3 pi 



240 



(4.68) 
(4.69) 



(4.70) 



f' 2 [54r 5 / (3) + 742r 4 /" - 175 



840 
1 

120 



//' [l26r 6 / (4) + 2226r 5 / (3) + 9170r 4 /" + 495 
r 2 ff" \8r 5 f i3) + 1376r 4 /" - 1935 



+— r 4 / 3 [r 3 / (4) - 20r 2 / (3) - 453rf ' + 696r 3 /' 2 - 1371/' + 99r 4 ff 
15 L 
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r 3 / 4 -25r 5 / (3) - 671r 4 /" + 58r 6 f 2 - 1741r 3 /' - 168 



+8r u f - 52r 9 / 7 rff" \Ur 5 f {3) - 33 



56 



When the radial function has the form f(r) 



„2a-2 



/(l +r ) (our Case I), it is possible 



to perform the radial integrals explicitly to get the result: 



1 + a 2 



6\a 



+vr(a 2 -4)(a 2 - 1 
Note that, for |a| 



lnm 



7T (24a 8 - 60a 6 + 11a 4 + 50a 2 - 25) 1 



1800a 6 sin(vr/|a|) 
90a 8 - 152a 6 + 553a 4 - 126a 2 - 280 1 



77T 



(4.71) 



11025a 8 sin(27r/|a|) m 4 
1 or 2, taking the limit |a| — > 1 or |a| — > 2 in the right hand side of 



(4.71) should be understood. (One may compare (4.71) with the corresponding form for 
the scalar effective action given in (3.8) of Ref. [2].) This large mass expansion result will 
be compared with the numerically determined effective action later. In Case II, it is not 
possible to obtain the associated radial integrals in a closed form, but we can evaluate them 
numerically. 



2. Numerically exact computation 

We now turn to our numerical evaluation method. First consider partial wave contribu- 



tions with I 7^ j (or I — j — 0). By solving the differential equations (4.2) numerically, we can 



determine the value for the ratio of two functional determinants according to the GY formula 



(4.1 ). One may easily solve the equation corresponding to the free equation. This free radial 
wave function is given in terms of the modified Bessel function, i.e., ^ (r) = l2i+i(mr)/r. 
As noted in Ref. [3J, it is convenient to consider the ratio of two functions 

KAr) = Ml (4-72) 



which has a finite value even though each of the numerator and the denominator diverges 
in the r — > oo limit. This ratio function TZij(r) satisfies the differential equation 

d 2 Tli ~ (I n (mr) \ dKt ,• , , m n tl „ n . 

—4f +[- + 2m^) — M - VijUtj = 0, (4.73) 

under the initial value boundary conditions 

n ld \ r=0 = l, 7^.| r=0 = 0. (4.74) 
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These differential equations share the same character as the ones one encounters in the 
evaluation of the scalar effective action studied extensively in Ref. [2]. 

In the present problem it is also necessary to evaluate functional determinants involving 
2x2 matrix differential operators, and for this we must solve the matrix differential equa- 



tions (4.7) with the boundary condition (4.6). Here again, instead of directly solving them 



numerically, we will consider a new matrix function 

7C;(r 



^f ree (r)' 

It satisfies the matrix differential equation of the form 

dlli 

dr 



dr 2 \r I 2 i+i{mr) J 



with the initial boundary conditions 



n. 



l\r=0 



1 
1 



l\r=0 



- VuKt = 0, 








(4.75) 



(4.76) 



(4.77) 



Then the functional determinant of matrix differential operator in ( |4.4 ) can be determined 
in terms of the ordinary determinant of the 2x2 matrix lZi(r = oo). 

Using the values TZij(r = oo) and TLi{r = oo) found by the above method, each group of 



the functional determinants in the right hand side of (2.38) can be numerically evaluated to 
find the value of T^^A; m): 



KL 



{A;m) 



(In ^o,o -ln72o,i) + J2 {(2/ + l) 2 (In det 7?., 

l=i,i,- 



(4.78) 



+ ln7£ ; _i i,i + hiTZ, , i ,_ i 

t 2'~ r 2 'i"2> 2 



ln^ M+ i -hiH l+ i 



Combining this with the contribution from the high partial wave part (keeping up to the 
terms of O(j^)) 



1>L 



(A;m)\ 



truncated 



dr 



Q 2 L 2 + Q X L + Q log In ( 2L ( U + 1 A + Q[ 



fir 



Q-i 0-2 

L L 2 



(4.79) 



we can evaluate the effective action to very high accuracy. In (4.79), explicit forms of 



Q 2 , ■ ■ ■ and Q can be found from ( |3. 18 ) and those of Q_i and Q_ 2 from Appendix |A| As 
explained at the end of Section |nj including the ^-suppressed terms Q-x/L and Q_ 2 /L 2 (of 
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m 



FIG. 4. Plots of the effective action in our Case I backgrounds, numerically evaluated as a function 
of m for various values of a (= 1, 2, 3, 4, 5). The solid lines (with the same color as that used for 
numerical values) are the corresponding results of large mass expansion. 



r[>i(^; m)) in our effective action formula (2.37) makes it possible to evaluate the effective 



action accurately with a relatively small value of L. In practice, with a choice of 20 < L < 50, 
we could obtain the value for the effective action with the accuracy of 10~ 6 . 

In Case I, the radial function is f(r) = r 2a ~ 2 / (l+r 2a ). In this case we have the expression 



r 



(3.17) for the high partial-wave part T l>L (A;m) and the Q-\ and Q-2 terms can also be 



evaluated. This high partial-wave part must be combined with the low partial-wave part 
r^^A; m) which requires extensive numerical work. We have evaluated the fermion effective 
action as a function of m for the values of a = 1,2,3,4,5 for concreteness. In Fig. El we 



plot these results for the full effective action r ren (74; m) (given by (2.31 )), together with the 



corresponding results based on the large mass expansion in (4.67). From this plot we see 



that the validity range of the large mass expansion varies with the 'stiffness' a (see Fig. 
[ija)) — it is valid if the mass m is such that mp > 1.2 for a = 1 (not stiff) and mp > 4 
for a = 5 (stiff). Note that our numerically determined curves for T ren (A; m) are essentially 
exact ones for all nonzero values of m. Further, one might even try to read the values of the 
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quantity linim^o [r ren + ln(mp)] from the extrapolation of our numerical data. (It is quite 
difficult to evaluate numerically the effective action with m = since some of the solutions 
to the GY equations diverge in the r — > oo limit). See TABLE [n] for the values of C(a) 
determined numerically for m = ^ (assuming the 'form' (4.58) against their exact values 



based on (4.59). The accuracy of the numerical method used by us is beyond question. 



a 


C(a) 


exact (m = 0) 


numerical(m = ^) 


1 


-0.291747 


-0.2916 


2 


-0.269189 


-0.2690 


3 


-0.378112 


-0.3782 


4 


-0.590437 


-0.5905 


5 


-0.883495 


-0.8835 



TABLE II. Exact values (with m = 0) versus numerical values (with m = 1/100) of C(a) for 
a = 1,2,3,4,5 



We now turn to Case II, where f(r) 



2 ( 1 ^ r 2) [1 + tanh(^^)]. There are two free pa- 



rameters (3 and R, with an instanton-like background (see Fig. pTc)) if /3 > and an 
instanton-antiinstanton configuration (see Fig. [T|(d)) if (3 < 0. For various choices of (3 
and R, we have determined numerically L ren (A;m), as a function of mass m. First, for a 
fixed value of R = 5 and (3 = 2, 1, | and |, the resulting functions are plotted in Fig. |5j 
Here notice that the instanton configuration becomes more stiff-wall-like if the value of (3 is 
reduced, and for the validity range of large mass expansion we find mp ^ 1 for (3 = 1 and 
mp > 4 for (3 = |. For the same R- value the fermion effective action L ren with negative 
values of (3(= —2, —1, — |, — |) are evaluated also, and they are plotted in Fig. [6j Clearly, 
as m approaches zero, the effective action becomes singular if > 0, but remains finite 
for (3 < 0; this is a phenomenon directly connected with the existence or nonexistence of a 
fermion zero mode in the m = system. 

We also studied how the fermion effective action changes as the parameter R is varied 
(taking here m = ^ and (3 = ±1): these results are in Fig. [7J From the plots shown in 
Fig. [7j one will notice that the fermion effective actions for both values of f3 = ±1 grow 
linearly with R if R becomes large. This should not be anything surprising — for the gauge 
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FIG. 5. Plots of r rcn (m) in our Case II backgrounds (with R = 5) for j3 = 2, 1, \ an d \- Near 
m = 0, all of blow up, exhibiting the ln(l/m) divergence. The solid lines (with the same color) 
correspond to the results of large mass expansion. 

background field involved here, the classical Yang-Mills action, which enters the effective 
action through renormalization counterterms, also grows linearly with R. (That the two 
curves in Fig. [7| one for /3 = 1 and the other appropriate to the case /3 = — 1, have even the 
same large- R slope is related to the point discussed below). 



Here recall that, according to the relation (2.22) and the remarks that follow immedi- 
ately, our Case II background with a negative value of (3 = —f3 (/3 > 0) can actually be 
viewed as a composite configuration involving an instanton located near the origin (which 
is gauge-equivalent to our Case I background with a = 1) and an antiinstanton-like config- 
uration associated to our Case II background with positive (3 = +/3q. Then, we may define 
the fermion-induced 'interaction energy' between the instanton and the (spherical-wall-like) 
antiinstanton, separated by distance R, as 



T int (R;m) = T^(R;m) 



Of=l 



m, 



(4.80) 



where the designation (I) or (II) refer to our Case I or Case II background, respectively. 
Choosing (3 = 1 for definiteness, we used our numerical results obtained for r ren (A;m) in 
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FIG. 6. Plots of r ren (m) for negative values of (3 (=—2,-1,-1/2, —1/4) when R = 5, together 
with the results of large mass expansion(represented by the solid lines). All of our numerical plots 
approach finite values as m — > 0. 



Case II background to study how this interaction energy depends on R, at some chosen 
values of fermion mass m. See our plots shown in Fig. [8] Clearly, with a finite value for m, 
this interaction energy vanishes as R becomes sufficiently large. But this interaction dies 
away at far slower rate as the mass value becomes very small — i.e., smaller the fermion 
mass, more long-ranged interaction seen between the instanton and antiinstanton. [As we 
remarked already, with strictly zero mass, it is very difficult to perform numerical study]. 

To see the origin of the above long-range interaction at small mass, it is useful to separate 
the contribution to the effective action r^(i?; m) (with /3 < 0) coming from the first brackets 



in the right hand side of (2.38) — the part containing the / = partial waves — from the 
rest. It is in this / = partial wave term (which we denote as T^ l=0 \R; m) below) where the 
effects related to the disappearance of normalizable fermion zero modes (at m = 0) in the 
instanton-antiinstanton composite configuration are relevant. Now see our numerical results 



for °^(R;m) (with /3 = —1) given in Fig. [9j and see also Fig. 10 where we present the 
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R 



FIG. 7. Plots of r ren for the mass value m = ^ as a function of R for two cases with /3 = 1 (blue 
dots) and ft = — 1 (red squares). 



plots for the effective action with the very / = contribution removed, i.e., for 

T int (R; m) = T int (R; m) - ^ l =°\R; m) 



(4.81) 



The fact that the latter quantity r int (i?;m) as a function of R becomes flat rapidly is an 
unambiguous sign that the I = partial wave term is mainly responsible for the above long 
range interaction between the instanton and the antiinstanton at small fermion mass. [The 
data for m = included here is the result of our direct calculation (using the exact massless 
GY wave functions) for the I = partial wave contribution, and at large R this curve is 



clearly consistent with the behavior found in (4.66) 



We found that the numerical data for the function °^(-R; m) presented in Fig. [9] are 
well approximated by the simple function 

' ,2 ! \ 



m 



ri'=° r L( J R;m) = -ln — + R , 



(A w 5.55). (4.82) 
Certainly, it is valid when m is small and R is large. As m — > it becomes the function 4 In R 



in (4.66), derived in Sec. IV A. In the R — > 00 limit, it approaches a constant, — him 2 /A. 
The values of - lnm 2 /A for m = 1/400, 1/200, 1/100, 1/50, and 1/10 are 13.69, 12.31, 
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m=l/10 
m=l/50 
m= 1/100 
m= 1/200 
m= 1/400 



FIG. 8. The interaction energy Ti nt as a function the distance R, for m = 1/10, 1/50, 1/100, 1/200, 
1/400. 
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10 20 30 40 50 60 

FIG. 9. Contribution from the I = partial wave term, ; m), as a function of R for various 

small values of m (=0, jgg, ^gg, gg, jg). The solid lines (with the appropriate color) represent 



our approximate formula in (4.82) for the given mass values 
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FIG. 10. The interaction energy with the very I = contribution removed, Ti nt (R;m). 
10.92, 9.53, and 6.31, respectively. These values, if the opposite sign is taken, are very 



close to the values denoted by flat lines in Fig. 10 (which correspond to the sum of all 
/ 7^ contributions to the interaction energy and the I = contributions to the effective 
actions for the separated instanton and antiinstanton configurations). Adding him 2 / A to 



the function in (4.82), we may thus get an approximate formula for the interaction energy 
in the form 



i(approx) 



int 



(R;m) = -lnfl + 



A 
w?R A 



(4.83) 



In Fig. yTJ plots for this function for the mass values chosen for our numerical works are 
given together with the related numerical data for comparison. Note that the numerical 
data and those our approximate formula coincide very well for R > 5, as long as m ^ 1/10. 



V. CONCLUDING REMARKS 

The partial-wave cutoff method has been extended to evaluate the 4-D spinor effective 
action in radially symmetric, non-Abelian, gauge backgrounds. With a suitable extension of 
our previous scheme, it is shown that this method retains its full power even when an infinite 
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FIG. 11. Plot of the same data in Fig. [8] together with our approximation function — ln(l + 
A/{m 2 R 4 )). 

number of functional determinants of matrix-valued radial differential operators have to be 
evaluated. By this method we have determined the fermion effective action numerically (for 
generic mass value) in various instanton-like and instanton-antiinstanton-like backgrounds. 
The validity range of large mass expansion has been checked using these numerically exact 
calculations. Also the results of our computation have been used to study the effective, light- 
fermion-induced, instanton-antiinstanton interaction. The technical aspects, elaborated in 
the present paper, should find useful application in other effective-action-related studies. 
One might also make use of our study as a basis to test how good non-Abelian derivative 
expansions [19] are. (Such study in the Abelian case was made in Ref. [5]). 

In the massless limit, evaluating the fermion effective action is in some sense simpler than 
evaluating the scalar effective action in the same background. This is due to the unique 
feature of factorizability that the massless fermion Gel'fand-Yaglom equations for partial 
waves have. By utilizing this factorization fully, we have in fact found analytically the small- 
mass-limit form of the full fermion effective action in our general Case I backgrounds. As one 
can always use the systematic large mass expansion result for the effective action if mass is 
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not so small, the additional knowledge on its small-mass-limit behavior is often sufficient to 
obtain a global fit to the full mass-dependence of the effective action [TH] . In this regard, some 
further effort might be desirable in trying to solve the massless fermion Gel'fand-Yaglom 
equations, with more general background functions than the forms considered here. 
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Appendix A: Subleading terms for the High partial- wave contribution T^^^Ajm) 



For the background field of the form ( l.lj ), (3.17) and (3.18) contain our calculated 
result for the high partial- wave contribution T^ L (A;m). But only the minimal terms were 
included in these formulas. In case we want to implement our effective action calculation 
scheme numerically, it will then be necessary to evaluate a very large number of partial-wave 
functional determinants from the low partial- wave side. But, if we include some -^-suppressed 
terms in our expression for the high partial-wave contribution, it becomes possible to take 
our partial wave cutoff L at relatively small value; that is, based on fewer calculations of 
functional determinants from the low partial- wave side, we can reach numerically convergent 
results for the full quantity v\y L {A\m). As we utilized this idea in our numerical analysis 
given in Sec. IV, we will here present the explicit forms of the 0( j-) and O(j^) terms which 



may be kept inside the integrand of our high partial-wave contribution formula (|3.17|). These 

-2( r )j2, are given by (here, u = 



terms, which we denote as Q_i{r)j- 
Q :ir) 



Q- 



Til mr 
1 ~f 1 2L 



32m 9 



6rV (21m 4 - 18m 2 + 5) f + 12rV (m 2 - l) //" - 6r V (u 2 - l) /" 
+12rV (27m 4 - 12m 2 + 5) f + r 2 (32m 8 - 389m 6 + 297m 4 - 75m 2 - 105) f 2 
-6r 4 M 4 (m 2 + l) if) 2 - 60rV (u 4 + l) //' - 2ru 2 (8m 6 + 3m 4 + 6m 2 - 15) /' 
+72r 5 M 6 / 2 /' + (-64m 8 + 191u 6 - 225m 4 - 15u 2 + 105) /] , (Al) 
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g_ 2 (r) = 768Qul l {u - 1} [-80m 4 {80m 6 , -(« + !) (103m 4 - 38m 2 + 7 ) } r 10 / 6 

+4800m 8 (2wm 2 - u - l) r 9 f 4 f - 240m 6 {4vm 4 + (u + 1) (llM 2 - 5) } r 8 / 2 (/') 2 

+240m 4 {i44m% - (u + 1) (99m 4 - 18m 2 + 7) } r 8 f 

+ 160m 6 {8m 4 V - (u + 1) (13m 2 - 5) } r 8 / 3 /" - 480m 1 V 7 {(/') 3 + / 2 / (3) } 

-160m 4 {ll2M% - 5(m - 1)(m + l) 2 (m 2 - 7)} r 7 ff - 1920u w vr 7 f f f" 

-240m 6 {36^m 4 + (m + 1) (7m 2 + 5) } r 6 / (ff + 24m 8 {2vu 2 + 3(m + 1)} r 6 (f'f 

-20m 2 {2496m 10 + (m + 1) (3576m 8 - 4377m 6 + 3135m 4 - 595m 2 - 315) } r 6 / 4 

-240m 6 {34m% - (m + 1) (m 2 - 5) } r 6 ff" + 48m 8 {3vu 2 + 2(m + 1)} r 6 f'f ( V 

+48m 8 (m + l)r 6 // (4) - 240m 4 {62m 8 + (m + 1) (8m 6 + 144m 4 - 5m 2 + 35) } r 5 ff 

-24m 8 (m + l)r 4 / (4) + 48m 6 {16wm 4 + (m + 1) (9m 2 + 10) } r 5 ff (3) 

+4m 4 {1208m 8 + (m + 1) (968m 6 + 707m 4 + 1190m 2 + 525) } r 4 {ff 

+40m 2 {384m 10 + (m + 1) (1800m 8 - 3561m 6 + 1935m 4 - 595m 2 - 315) } r 4 f 

+8m 4 {864m 8 + (m + 1) (1104m 6 - 489m 4 + 900m 2 + 385)} r A ff" 

+ 16m 6 {139wm 4 + (m + 1) (76m 2 + 55) } r 5 f'f" - 40m 6 (m + 1) (2m 4 + 3m 2 + 6) r 3 / (3) 

+40m 2 {288m 10 + (m + 1) (24m 8 + 657m 6 - 597m 4 + 973m 2 + 315) } r 3 ff 

-15(m + 1) (3072m 10 - 7779m 8 + 3820m 6 + 4270m 4 - 4788m 2 - 1155) r 2 f 2 

-20m 4 (m + 1) (64m 6 - 159m 4 + 126m 2 + 77) r 2 f" 

+60m 2 (m + 1) (14m 8 - 41m 6 + 327m 4 - 231m 2 - 105) rf 

+15(m + 1) (1040m 10 - 2823m 8 - 400m 6 + 6930m 4 - 3528m 2 - 1155) /] , (A2) 

where v = m 2 + u + 1 and /', /" and denote the first, second and n-th derivatives of 
f(r). 

Appendix B: Small mass limits of Gel'fand-Yaglom wave functions 

In this appendix we will describe the method that allows to determine the asymptotics of 
massive GY solutions when the mass is sufficiently small. We will here assume that the exact 
massless GY solution is known. Then, depending on whether the massless GY solution is 
normalizable or not, we have two different ways to obtain the desired GY solution for small 
nonzero mass with global validity (i.e., the asymptotic region included). 
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First, suppose the exact m = solution does not correspond to a normalizable zero mode 



(like the GY solutions given in (4.22), (4.23) and (4.44)). Then we can follow the idea of 
Ref. [32] and write the perturbative GY solution of (4.2) for small nonzero mass as 



(Bl) 



where ipo{ r ) is the known massless solution, and ipi(r), ^(Oj ' ' ' denote appropriate mass- 
independent functions. We will refer to this solution as the small-r solution. Being concerned 



with the leading small mass behavior, we will here keep the leading order solution in (Bl) 



only. This naive solution works fine when the mass m is very small and at the same time 
mr can be taken to be finite. But, in our case, we are interested in the large-r asymptotic 
behavior of i/)(r;m), i.e., at r ^> — . Therefore we need to consider another perturbative 
solution which we call as the the large-r solution. For this large-r solution we change the 
variable r to x = mr (see Ref. [IB] for more detailed discussions) and write the corresponding 



solution as (p(x) (instead of ip(r)). We can then recast the GY equation (4.2) as 

r d 2 



dx 2 



3 d 41(1 + 1) 1 , fx 

-5-+ 2 - + -2^ ~ 
x ox x z m z \m 



0. 



(B2) 



From (2.26)-(2.29) and assuming that the function H(r) = r f(r) approaches 1 as r — > 00 



(i.e., for Case I and also for Case II with (3 > 0), we can approximate the potential in (B2) 
such that 



4Z(Z + I) 



x A 



1 

m 2 



Aq(q + 1) 



0(m' 



(B3) 



(q denotes the quantum number defined in (2.23)). Note that this approximation is good as 



long as — is not so small(, say, compared to 1). With only the first term in (B3) kept, we 



have the leading order solution of (B2) in the form 



M%) = (const.)2(2g + 1)! f-V" 

\mj x 



(B4) 



which is proportional to the free solution of (B2) with / replaced by q. The proportionality 



constant will be determined at some intermediate point r = R (or, equivalently, at x — mR), 
for R satisfying the condition 1 R — . The small-r solution should be valid for 
< r < R, while the large-r solution is valid for r > R (i.e., x > mR); hence, both solutions 
are valid near the point r ~ R, so that we can demand 



ip(R) = <p{mR). 



(B5) 
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If if)o(R) has an asymptotic behavior ^q(R) ~ cqR 2q (with an appropriate constant cq), 
we can now fix the overall constant in (B4) (i.e., demand (po(mR) ~ CoR 2q ) to obtain the 
complete, zeroth order, solution 

M^) = 2c (2 q +1)\(-) 29I ^1. (B6) 
\mj x 



Based on the form (B6), we can obtain the correct asymptotic behavior of the GY wave 
function: 



l~2 ( 2 \ 2q e x 
Hm if)(r,m) ~ lim (p (x) ~ cm/- — {2q + l)\—j-. (B7) 

1 — >oo 1^00 \ 7T \TTl / X I 



Consequently, for the ratio of GY solutions, we find 



— ^fe»(r;m) (2Z + 1)! VW v ; 



T— >OQ 



Note that, as m goes to zero, the above ratio goes to infinity or zero or remain finite 
depending on the value of q. Let us apply this to our Case I. The functional determinants 
for partial waves corresponding to j — I ± 1 at small mass limit can be found immediately. 



From (4.22) and (4.23), the large- R behaviors of the massless solutions are 



VWr;m = 0)~i? im , (B9) 
21 + 1 

*(; l>l _ 1 (r;m = 0)~^ r R 2l -\ (BIO) 



Thus, by identifying cq from these behaviors and plugging them to (B8) (with q = I ± | for 
j — I ± 1), we find that 

det(-H M+1 +m 2 ) ^ M +i(r;m) 4(1 + 1) 

— — -—r ^— = lim —? — ~ , (Bll) 

det("Hf ec + m 2 ) r ->°° %jjf ee (r; m) m 

det(Hi,i-i + m 2 ) ipu-i(r;m) m , . 

— 1 = hm — i — f ~ — • B12 

det(^[ ree + m 2 ) r^oo ^f ree (r;m) 41 y ' 

Now one can verify explicitly that these small mass limits yield, if used for a certain specific 
group (indicated in the main text), the same results as those based on the exact zero mass 



results such as (4.27) with the same combination. 

For j = I = (and when the function H(r) = r 2 f(r) approaches 1 as r — > 00), it is 
rather difficult to follow the method of Ref. |18|. In this case, the massless GY solution 



(4.29) corresponds to the normalizable zero mode. As a result, the small-r solution cannot be 



matched with the large-r solution of the form (B4). Here we need to consider the Macdonald 



function for the large-r solution; but, the Macdonald function vanishes exponentially as r — > 
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oo, and therefore we need to take into account the next order term which is mathematically 
very complicated. In this situation the method of Ref. [I], which utilizes the normalizable 
character of the massless solution in a crucial way, can be simpler. Note that from the 
massive and massless GY equations 

id 2 3d 4Z(Z + 1) 



+ 2 + V<j(r) + mf \ V>(r; m) = 0, (B13) 



° 2 39 + 4 ^ + V lj (r)U, (^0, (B14) 



we can deduce the relation 

d 

— \r 3 ?p'(r; m)ipo(r) — r 3 ip(r; mjipLir) \ = m 2 r 3 ip(r; m)^o(V), (B15) 
or L > 

Integrating this equation from r = to r = R e and then dividing by Rjip' (R e ), we get 

^ R - m) few - rJ = m — mm — • (B16) 

For large R e , we can write ip Q (R e ) ~ ^ (as H(r) approaches 1 as r — > oo) and ip(R;m) ~ 
Be mRe , where A and B are some constants. Hence, by considering the R e — > oo limit with 



(B16), we are led to conclude that 



fn Re r 3 ib(r:m)ibn(r)dr ,„ 
lim ij(R e ;m) = -3m lim ^ ^ • B1 ? 



When mass m is small, this equation simplifies to 



lim ^(R e ; m) ~ -3m , ii^^M^- , (B18) 



For our Case I the explicit massless solution is available from (4.29): 



^ ,o(r; m = 0) = Mr) = (B19) 

_|_ ]\ 



Using this solution, we can evaluate the expression appearing in the right hand side of (B18) 

explicitly. The result is finite, and in this way we can secure the following result for the 

related functional determinant: 

det(H ,o + m 2 ) = ^(R e ;m) ^ r(l + |) r(|) 

det(-<9L + m 2 ) R^oo^™(R e ;m) 2T( 3 ^ 



'(0) 

We now turn to the case with j = I ^ 0, where the GY solutions are given by 2 x 2 



matrices. First, note that we can approximate the potential in (B2) for small mass m as 



x z m z \mj x z I q 
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Thus the leading order solution of (B2) using the matrix potential given above is 

/ 



2(2/ + 2)! 



21+1 







V 



o 



2(20! (|) 



21-1 



Co, 



(B22) 



where Co is a constant 2x2 matrix which can be determined using the condition (B5) at 
r = R. To find Co, we need to know an asymptotic behavior of the massless GY solution. For 



this j = I 7^ the massless solution is given in (4.44), with relevant expressions also in (4.45), 



(4.42) and (4.46). Since this solution was found after a suitable unitary transformation was 



performed (with the help of the unitary matrix U in (4.33)), our condition (B5) translates 
to 

' R 2l+1 



U^i(R)U ~ <p (mR) 



R 



21-1 



Cq. 



(B23) 



As the asymptotic behavior of the determinant of *&i(R) is given in (4.52), we find that 



det c 



o 



det*,(i2) (2/ + l) 3 r(^ 



, 2 ' 



(B24) 



R41 8j(i + i) 2 r(f)'r(^y 

The true asymptotic behavior of the determinant of the GY wave function can then be found 



as 



lim det if) u(r;m) ~ lim det<^ (x) 



2 e zx f 2 



ii 



it x 3 \m 

with det c given by ( B24[ ). Based on this, we obtain the following result: 
det(^ M + m 2 ) detV M (r;m) 21 + 2 

detc . 



(2/ + 2)!(2/)!detc (B25) 



detV M (r;m) 
lim ■ 

r-*x> ijjf ee (r;m) 2 



(B26) 



{det(-9 ( 2 /)+ m2)} Z ™ *Pr(r;my 21 + 1 
Here note that det Cq corresponds to the result of the related massless functional determinant. 



Therefore (B26) apparently shows that the small mass limit of the functional determinant 



is different from the massless functional determinant. Only with the combination in (2.38), 
the small mass limit and exact zero mass analysis yield the same results. 
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